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Abstract 

In ultrarelativistic heavy ion collisions a finite size QCD medium is created. In this paper we 
compute radiative energy loss to zeroth order in opacity by taking into account finite size effects. 
Transition radiation occurs on the boundary between the finite size medium and the vacuum, and 
we show that it lowers the difference between medium and vacuum zeroth order radiative energy 
loss relative to the infinite size medium case. Further, in all previous computations of light parton 
radiation to zeroth order in opacity, there was a divergence caused by the fact that the energy loss 
is infinite in the vacuum and finite in the QCD medium. We show that this infinite discontinuity 
is naturally regulated by including the transition radiation. 


1 Introduction 

The suppression pattern of high transverse momentum hadrons is considered to be a powerful tool to 
map out the density of the produced QCD plasma m-M- This suppression (so called jet quenching) 
is assumed to be mainly due to the medium induced radiative energy loss of high energy partons 
propagating through ultra-dense QCD matter ail!- However, even if final state multiple elastic and 
inelastic interactions are neglected, the difference between medium and the vacuum 0*^ order energy 
loss would still be significant. This is due to the fact that the gluon dispersion relation is different 
in the medium and the vacuum, leading to differences in the associated 0*^ order radiation. This 
effect was first pointed out by Ter-Mikayelian ElEI, who considered the QED plasma case. In dniiiii 
we developed a non-abelian QCD analog of the Ter-Mikayelian plasmon effect for the case of heavy 
quarks. We showed that while the Ter-Mikayelian effect is negligible for bottom quarks, it has an 
important effect on charm quarks, since it leads to a significant reduction of the vacuum radiation. 
This result is a consequence of the fact that the gluons in the QCD medium acquire a finite mass 
proportional to the temperature of the medium. 

The computation presented in cni was done under the assumption of an infinite QCD medium. 
This raises the question how to generalize these results to the more realistic case of a finite size QCD 
medium created in ultra-relativistic heavy ion collisions (URHIC). How the results from will be 
modified by finite size QCD effects is the first goal of this paper. 

Further, it is well known that light parton 0*^ order energy loss is not infrared safe, i.e. it goes 
to infinity when the parton mass goes to zero. This infrared divergence is absorbed in the DGLAP 
evolution m, so that the only part which contributes to the jet quenching is the difference between 
medium and vacuum energy loss. Since all previous computations [ilZl assumed that the light quarks 
and gluons have the same zero mass in both the medium and the vacuum, the difference between 
medium and vacuum energy loss was found to be finite. 
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However, with the introduction of the Ter-Mikayelian effect, the finite parton mass in the medium 
regulates the infrared divergence of the 0*^ order energy loss in the medium, while the corresponding 
energy loss in the vacuum remains infinite. This leads to the question how to regulate this discontinuity 
between medium and vacuum light parton energy losses. The second goal of this paper is to show how 
transition radiation naturally solves this problem. 

The outline of the paper is as follows. In Section |21 we will compute the 0*^ order radiation in a 
finite size QCD medium for both light and heavy quarks. For charm quarks we will show that the 
transition radiation lowers the Ter-Mikayelian effect from 30% m to 15 — 20%. Additionally, we 
will show that for light partons the transition radiation naturally regulates the infinite discontinuity 
between 0*^ order medium and the vacuum radiative energy loss. In Section |31 we will study how 
the difference between medium and the vacuum 0*^ order energy loss depends on where the particle 
is produced. We will show that the difference between medium and vacuum 0*^ order energy loss is 
positive (as intuitively expected) as long as the probe is produced far outside the medium (QED case). 
However, if the particle is produced inside the medium, such as in the QCD case, we will obtain that 
this naive expectation may not hold. In Section we will extend our study from Sections |21 and |21 
to include the fact that due to the confinement in the vacuum, the gluons may acquire finite mass. 
We will obtain qualitatively different results, compared to those presented in Sections and 01 if the 
gluon mass in the vacuum is larger than in the medium. In Section 0 we will combine the results 
presented in Sections [21 and |11 with the medium induced radiative energy loss |13| . We find that for 
certain realistic values of the gluon mass in the vacuum, the light quarks can leave the L < 3 fm 
medium essentially unquenched. We argue that the results presented in the Section 01 may provide us 
with a hint toward solving the puzzle posed by M- Finally, in Section 01 we summarize our results 
and put our work in the context of future research. 

2 The one gluon 0^^ order radiation in a finite size QCD medium 

The aim of this section is to compute the 0*^ order radiative energy loss when the parton is produced 
in a finite size dielectric medium. To introduce the finite size medium, we start from the approach 
described in m- As in HSl, we consider the static medium of size L, and define two gluon masses, 
(for gluon radiated in the vacuum), and rug^p (for gluon radiated in the medium). We also assume 
that, in general, the running coupling constant can be different in the vacuum and in QGP. However, 
contrary to m, we ignore spin effects, since they are irrelevant in the soft radiation limit that we 
consider in this paper. 

To compute the 0*^ order radiative energy loss in a finite QCD medium, we have to compute the 
squared amplitude of a Feynman diagram, M™'^. The Feynman diagram represents the source J, 
which at time xq produces an off-shell jet with momentum p' and subsequently (at xi > xq) radiates 
a gluon with momentum k. The jet emerges with momentum p and mass M. We neglect the thermal 
shifts of the quark mass since 1) for heavy quarks, thermal effects on the quark mass are negligible, 
and 2) light quarks will be treated as massless particles for the reason explained in footnote 2. 

The matrix element for this 0*^ order in opacity radiation process can then be written in the 
following form 

j^rad ^ J ^4^ ^4^^ Am(xi - Xq) u'^(xi) Aj,(xi)^>'^(xi) (1) 

where <h(xi) = is the wave function of the final quark with (on-shell) momentum p and A^(xi) is 
the wave function of the emitted gluon. Vertex function u^(xi) is given by u^(xi) = ( 7 (xi)( — d^), 

where ^(xi) is the running coupling constant which is in general different in the medium than in the 
vacuum. 
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For this problem it is convenient to use light cone coordinates m- This coordinate system is 
appropriate for systems moving with almost the speed of light. It is obtained by choosing new space- 
time coordinates [x"*",x“,x],^ related to the coordinates in the laboratory frame (t,z,x) by (x is the 
transverse coordinate) 

x^ = {t + z), x~ = {t — z). (2) 

In the same way the light cone momentum [p'^,p~,p] is related to the momentum in the laboratory 
frame [E,pz,p] by (p is the transverse momentum) 

P'^ = {E + Pz), p~ = {E-p^). (3) 

Additionally, it can be shown that in the light cone coordinate system the propagator Am{x) reduces 
to (see [in]): 


=(2^ / (4) 

where p'~ = ^ . 

In the spinless case, the wave function of the emitted gluon with momentum k can be written as 

A^{x) = ei,{k) ^>g(x) c, (5) 

where e(fc) = [o, 2 |^,e] is the transverse polarization and c is the color factor of the radiated gluon. 


^>g(x) = e 


—i|[fc+x +f d^k (^)]+ik-; 


( 6 ) 


is the wave function (derived in Appendix A) that satisfies the Klein-Gordon equation with position 
dependent gluon mass and A:“(x^) =- ^ -. In the static approximation m^(x+) becomes 


X 


x^ 


m, 


s(x ) = mg^p e{L - ^) + 'mg,v 0{-Y - 


( 7 ) 


where rug^p is gluon mass for the gluon radiated in the medium, while rug^y is gluon mass for the gluon 
radiated in the vacuum. 

We can now compute M^ad by substituting Eqs. 0-0 in Eq. o, which leads to the fallowing 
result (see Appendix B for detailed calculation of M^ad)'- 


€*k 

Mrad = -2iJ{p + k) - 

X 


Qp 1 - Qy 


p. 


/+ 




Xv 


( 8 ) 


where Pp [py) is the running coupling constant in the medium (vacuum). The variable x is defined as 


= and 

n'-r 


Xv = 

Xp = 


+ TTig „ 

xp'+ 


xp 


/+ 


(9) 


^Note that the x"*" and x axes of the new frame he on the light cone m- 


3 













Here, we use for the initial quark a plane wave state in the x-plane and set = 0. Then p = —k, 

_ /\/f2 , v2_LA/f2 


^ and p = c^_'tw+ ■ Note that soft radiation is defined as x <C 1 (i.e. p'^ 3> k~^), so we 


1? Additionally, as in [m, we assume that J varies slowly with p, so that 


assume that 1 — x 
J{p + k) K J{p). 

In soft gluon approximation, the spectrum can be extracted from Eq 


as (see jTT]) 


i2 

2E(27r)3 2cu(27r)3 


d^Njd^N^ , 


where 


d^Nj = dR\Jip)\^ 


d^p 

(27r)32p‘^ 


( 10 ) 


( 11 ) 


Here dji = 3 (for three dimensional representation of the quarks). 

Finally, by nsing Eqs. (|Tni[THl together with Eq. (jH)) we obtain the main order fractional energy 
loss (/ = AE/E) for massive quarks and gluons in the QCD medium of finite size L 


dl 
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dx d?'k 
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dl^vl 

dx 
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+ 2 


Cnk^ V«s(^) ( 


7r2(p'+)2 

CrOLS 


Xp 


Xp 


Xv 


(1 - cosixpL)) 


- ml p) 


dx d2k 7r2 (k2 + m2p + M2x2)2(k2 + m2^^ + M2x2) 

(k^ + m^ +m2x2)L 

x(l — cos(- — -)), 

^ ^ 2Ex 


( 12 ) 


where C(^{k) {ct'g{k)) is the running coupling constant in the medium (vacuum). Cr is the color 
Casimir for the partons, i.e. Cr = 4/3 for quarks and Cr = 3 for gluons. E is the initial jet energy. 
The second equation in m is valid in the o?g{k) = a'g{k) case, where the physical meaning of the 
obtained results are more evident. 

is 0*^ order in opacity fractional energy loss in the vacuum: 


dlil CRa^jk) k2 

dx d2k 7r2 (k2 + m'^ ,^ + M'^x'^Y 


and is total 0*^ 


order in opacity fractional energy loss in finite size medium, given by 


r( 0 ) r( 0 ) I r( 0 ) 

^med ^TM ' ^trans ■ 

Here, Ir'Ij is the 0*^ order in opacity fractional energy loss in the infinite size medium, which can be 

obtained from by replacing mg^v by rug^p in Eq. m- ^trans the additional transition radiation 
occurring when the jet is traversing from the medium to the vacuum. As a crosscheck, we note that, 
by neglecting spin effects in ^3, Eq. (9) from can be reduced to the Eq. (IT^ above. We note that 
the computation in US] was done in 3-dimensional coordinate space {z, x), while our computations 
were more consistently performed in the light-cone 4-dimensional coordinate space. 

^Note that in the light quark case we set M = 0 GeV. Keeping the finite fight quark mass in Eq. 0 would correspond 
to retaining small corrections, which would be inconsistent with 1 — x ~ 1. 
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Figure 1: The reduction of the fractional 0*^ order energy loss for charm quarks due to both the 
QCD Ter-Mikayelian effect and transition radiation is shown as a function of the quark energy. The 
dashed-dotted curve shows the vacuum energy loss if gluons are treated as massless and transversely 
polarized. Dashed curve shows the effect using the Ter-Mikayelian effect only, and solid curve shows the 
total effect by including both the Ter-Mikayelian effect and transition radiation. Assumed thickness 
of the medium is L = 5 fm, charm mass is M = 1.5 GeV and the gluon mass is nig^v = 0 GeV 
{rrig^p = 0.35 GeV) for the vacuum (medium) case. 

To obtain the fractional energy losses, we perform the integration by using 0 < |k| < 2x{\ — x)E, 
i.e. in our computations we do not introduce a lower momentum cutoff. For running coupling we use 
the ’’Frozen a model” m 

47T 47T 

-M^) = Min{0.5, ——= Min{0.5,- 1,2 ,^2 ,^^ 2^2 } , (15) 

QCD 

where /^o = ^ for effective number of flavors n/ ^ 2.5, Aqcd ~ 0.2 GeV, and — M"^ = ^ ^ ^ 

Note that Og (a^) is obtained by setting = rrig p {rrig = „) in Eq. m- 

We next use Eqs. (|T^-([T^ to provide a comprehensive analysis of the influence of the transition 
radiation to the total 0*^ order in opacity energy loss for both light and heavy quarks. We first 
concentrate on the charm quark case and look how the 0*^ order energy loss depend on the initial jet 
energy (Fig.^ and thickness of the medium (Fig. El) • In Fig.^we see that for charm quarks, transition 
radiation lowers the Ter-Mikayelian effect from 30% to 15 — 20% for L = 5 fm medium. In Fig. Ejwe 
see that for a medium thickness greater than 4 fm the transition radiation becomes approximately 
independent of the thickness of the medium. 

The previous two figures were computed by assuming running coupling (given by Eq. (I15j) l which 
is different in the medium and in the vacuum. In Eig. Elwe want to test how the obtained medium and 
the vacuum fractional energy loss difference — Imc) is robust against variations in the choice of 

coupling constant. To do that, we plot as a function of initial jet energy for three different 

choices of running coupling as well as constant coupling as = 0.3. We see that for heavy (c and 
b) quarks, the difference between medium and the vacuum energy loss is almost independent on the 
choice of coupling constant. Eor the light quark case, we see that while is robust to the 

choice of running coupling, it is fairly sensitive to the choice between running and constant coupling. 


5 









0.35 


0.3 


AE 


( 0 ) 

vac 


o 

Hi 

<] 


0.25 


0.2 


A 


0 1 2 3 4 5 6 7 

L[fm] 

Figure 2; The reduction of the fractional 0*^ order energy loss for 10 GeV charm jet due to both 
the QCD Ter-Mikayelian effect and transition radiation is shown as a function of the thickness of the 
medium. The dashed-dotted curve shows the vacuum energy loss if gluons are treated as massless 
and transversely polarized. Dashed curve shows the effect using the Ter-Mikayelian effect only, and 
solid curve shows the total effect by including both the Ter-Mikayelian effect and transition radiation. 
We assume that the mass of the charm jet is M = 1.5 GeV and the gluon mass is rng^v = 0 GeV 
= 0.35 GeV) for the vacuum (medium) case. 




Figure 3: The difference between medium and the vacuum 0*^ order fractional energy loss for light 
(the lower four curves), charm (the middle four curves) and bottom quarks (the upper four curves) 
is shown as a function of initial jet energy. Solid curves are computed by assuming running coupling 
given by Eq. (USD, with different ofg in the medium and Og in the vacuum. Dashed, dot-dashed and 
dotted curves are computed by assuming the same coupling constant both in the medium and in the 
vacuum. For dashed (dot-dashed) curves we used running coupling oPg {ocg) given by Eq. (1151) . Dotted 
curves are computed by assuming constant coupling as = 0.3. Assumed thickness of the medium is 
L = 5 fm, light quark mass is M = 0 GeV, charm mass M = 1.5 GeV, and bottom mass M = 4.5 GeV. 
The gluon mass is = 0 GeV {rrig^p = 0.35 GeV) for the vacuum (medium) case. 
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Figure 4; The difference between medium and the vacuum 0*^ order fractional energy loss for light 
(dotted curve), charm (dashed curve) and bottom quark (dot-dashed curve) is shown as a function 
of the thickness of the medium. The initial jet energy is i? = 15 GeV. The curves were computed 
by assuming running coupling given by Eq. Hi, with different coupling in the medium and in the 
vacuum. Light quark mass is M = 0 GeV, charm mass M = 1.5 GeV, and bottom mass M = 4.5 GeV. 
The gluon mass is = 0 GeV {jUg^p = 0.35 GeV) for the vacuum (medium) case. 

For example, in p_\_ < 10 GeV range, 30% smaller — I^dc is obtained when constant coupling 
as = 0.3 is employed. 

Additionally, from Fig. |31we see that the finite mass (a.k.a. dead cone ^Hl) effect on 4^]^ — I^ac 
is strong, i.e. we see a qualitative difference between light and heavy quark 0*^ energy losses, which 
persists at high momentum. For example, while 4®^ — Ivac for light quarks is large and shows a 
noticeable dependence on jet energy, it is negligible for bottom quarks in the whole jet energy range. 

The most striking observation from Fig. |21is that, for the light quark case, the difference between 
medium and the vacuum energy loss (4^g^—4^ac) is finite. To validate this numerical result analytically, 
we assume a perturbative vacuum (i.e. rug^^ = 0 GeV) and the same coupling in the medium and in 
the vacuum (i.e. = a'g). In the (M = 0 GeV) light quark case, the Eq. (IT^ reduces to 



d(4, 


( 0 ) 

med 


— 4 ^°^) 

^vac) 


dx d^'k 


= -2 


Cnas 


m 


9,P 


vr^ (k2 -h mj pY^ 


1 — cos( 


(k^ -b 


9,P 




2Ex 


(16) 


This equation is infrared safe when k — > 0. This is an important result having in mind that without 
transition radiation, the Ter-Mikayelian effect leads to a discontinuity between finite medium and infi¬ 
nite vacuum energy loss. Therefore, we conclude that the transition radiation has special importance 
in the case of the light quarks since it provides a natural regularization of medium dispersion effects. 

In Fig. |l]we fix the energy jet to 15 GeV, and test how the idUd ~ depends on L for different 
types of quarks. We see that light and heavy quarks show significantly different thickness dependence. 
For heavy quarks, 4^]^—4.^ac saturate after some value of L (i.e. after L = 4 fm for charm and L = 1 fm 
for bottom). This saturation behavior is expected having in mind that for massive quarks and large 
enough L, cos(xpL) becomes rapidly oscillating function. Then (cos(xpL)) —> 0, and Eq. (IT^ becomes 
independent on the thickness of the medium. For the light quarks, we see that |4^g^ — I^acl increases 
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approximately linearly with L. This linear thickness dependence persists for higher jet energies as 
well (results now shown). This result is unexpected, having in mind that, in the light quark case 
and L/E ^ 0 limit, by expanding cosine in Eq. (jlbl) . we expect to obtain quadratic {Lp‘) thickness 
dependence for — Imc- However, by integrating Eq. (tTKl) in the L/E ^ limit (and without 
expanding cosine) we obtain a result which depends linearly on the thickness L of the medium (in 
agreement with the numerical results shown in Fig. El): 


/(O) _/(0) 

med ^vac 


CRasKp^ ^. E 1 

2 E ‘^^g,p 


(17) 


Finally, from Figs.^Elwe see that the total energy loss in the medium is smaller than in the vacuum. 
This result comes from Eq. (d, where we see that AE^^ — AEy^^ oc (ur^ ,, — rrig^p)- Therefore, if 
the gluon mass in the medium is larger than in the vacuum, then the total 0*^ order energy loss in 
the medium will be smaller than in the vacuum. Though mathematically correct, this result seems 
surprising, since it would be expected that the radiation in the medium (even at the lowest order) is 
always larger than the corresponding radiation in the vacuum, particularly having in mind the work 
presented in 0 Eni- In 0 EHI the transition radiation was studied for the particle traversing the 
QED medium, and it was shown that the lowest order radiative energy loss in the medium is always 
larger than in the vacuum. This study considered the case when the particle is produced outside the 
medium (at zq = —oo). However, contrary to the usual experiments involving a QED medium, where 
the medium is probed by using test particles produced far outside the medium, in URHIC the probes 
are produced inside the medium. Therefore, to be able to intuitively understand the result obtained 
from Eq. m, we have to take into account the qualitative change in the experimental approaches 
between QED and QCD. In the next Section we study how the difference between medium and the 
vacuum energy loss depends on the point of particle production. 
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Dependence of the 0^^ order energy loss on the point of particle 
production 



Figure 5: An illustration of the system studied in this Section. Gray area shows the medium of 
thickness L. Probe is produced in the vacuum at distance |/o| from the medium. 

To study this particular problem we assume the static system shown in Fig. El That is, we consider 
the probe produced in the vacuum at the finite distance Iq (i.e. x(J" = —21q < 0) from the medium of 
size L. As in the previous section, we assume that a gluon is subsequently radiated at a point xi. A 












gluon radiated in the vacuum (medium) has the mass rug^y (mg^p). Therefore, the gluon mass nig^K^) 
has the following form 


x~ x~ 

mg{-xl) = rug^y 6 ^(“) + rug^p e{-x.f) e{L - -^) + mg^y 9{^ - L). 


(18) 


As in the Section |21 we can now compute Myad-, by substituting Eqs. 0-® and m in Eq.(| 
which leads to the fallowing result (for the derivation see Appendix C): 


ek 1 

9v 

9v 

9p 
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.Xv 

Xv 

Xp 


Mrad = -2iJ{p + k)- 
Eq. (HUll together with the Eqs. (HTTl) and (tTTl) . leads to 


(19) 


^y^med 


dx 


_Cb^ ( \I^W 

7r2(y+)2 ^ 
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(1 - COS(xpL)) 


^2(y+)2 ^ 



(cos(y^Zo) 


cos(y^Zo + yp^))-(20) 


Equation CT represents the difference between medium and vacuum fractional energy loss when 
the probe is produced outside the medium. It is useful to look at two important limits of this equation: 

1) Iq —> 0. In this limit we recover the case when the particle is produced in the medium of size L. 
In this case, the Eq. CT reduces to Eq. m from the previous section. 

2) Iq —> oo. This limit corresponds to the case when the particle is produced far outside the 
medium, i.e. it is equivalent to the QED case studied in [^. 

When Iq ^ oo the second term in Eq. CT goes to zero. For ag{k) = ag{k) we obtain a result 
which agrees with [Hj: 


^y-^raed | 

dx d^k 


yP —(yV 

^S~^S 


,CRas{k) 


k^( 


- fnl,yf 


vr^ (k^ + + m?gyY{k^ + M^x^ + 'm’^p)'^ 


X [1 — cos( 


(k^ + M^x^ + mgp)L^ 


2Ex 


( 21 ) 


We see that, when Iq ^ oo, Eq. CT becomes positive definite independently of the gluon masses. 
Therefore, in this case the 0*^ order energy loss in the medium is always larger than in the vacuum. 
This result agrees with our intuitive expectations and with [^j. 


Finally, in Fig. IHlwe show the difference between medium and vacuum 0*^ order energy loss as a 
function of coordinate of particle production. The figure clearly shows the transition from positive 
definite values (for the case when the particle is produced far outside the medium) to negative values 
(obtained when the particle is produced inside the medium). 
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Figure 6 : The difference between medium and the vacuum 0*^ order energy loss for charm quarks 
is shown as a function of coordinate of particle production. The curves are computed by assuming 
running coupling given by Eq. Hi, with different couplings in the medium and in the vacuum. The 
thickness of the medium is L = 5 fm. The charm mass is M = 1.5 GeV and the gluon mass is 
'nig,v = 0 GeV {mg^p = 0.35 GeV) for the vacuum (medium) case. 

4 Dependence of transition radiation on the vacuum gluon mass 


The analysis presented in the previous two sections is based on the assumption that rug^^ = 0 in 
URHIG. However, this is true only in the case of the perturbative vacuum which does not take 
confinement into account. A phenomenological way to simulate confinement in the vacuum is to 
introduce an effective gluon mass nig^v 7 ^ 0 m In this case, there are two different vacuum gluon 
masses that can be found in literature |22j-|21], i.e. rug^y w ^qcd and rug^y ^ 0.7 GeV. 

To test how the difference between medium and vacuum 0*^ order energy loss depends on the 
different choices of rug^y we first show alternatives of Figs. Eia andinifor the case of rrig^y ^ ^qcd and 
rug^y ^ 0.7 GeV (see Figs. EUHl and El respectively). The left panels of Figs. |3 |H1 and El correspond to 
rug^y K 3 A.QCD case. We see that these figures are qualitatively similar to the rug^y = 0 case, although 
the net effect is smaller for nig^y « Aqcd- This result is expected, since oc (m^^p — m?g y), 

and Aqcd > 0 GeV. On the other hand, the results shown in the right panels of Figs. Q El and El are 
qualitatively different from the rug^y = 0 case. This is due to the fact that, in this case, the gluon 
mass in the vacuum is larger than in the medium (nig^y = 0.7 GeV), and thus we would expect that 
— ly^'^ would always be positive definite, in agreement with these figures. 

In Fig. uni we show the difference between medium and the vacuum 0 *^ order energy loss for light, 
charm and bottom quarks as a function of the vacuum gluon mass {nig^y). From this figure we see that 
for bottom quarks the difference is negligible. For charm quarks, in the range of experimental interest 
(0 GeV < mg^y < 0.7 GeV) this difference is in the range of about ±5%. In ESI EH] we showed that 
this difference has small effects on the heavy flavor experimental observables (i.e. negligible effect on 
bottom, and less than ± 0.1 on charm Raa)- However, for the light quarks, we see that the difference 
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Figure 7: The difference between medium and the vacuum 0*^ order fractional energy loss for light 
(dotted curve), charm (dashed curve) and bottom quark (dot-dashed curve) is shown as a function 
of the jet energy. The thickness of the medium is L = 5 fm. The curves are computed by assuming 
running coupling given by Eq. (USD, with different couplings in the medium and in the vacuum. Light 
quark mass is M = 0 GeV, charm mass M = 1.5 GeV, and bottom mass M = 4.5 GeV. The gluon 
mass in the medium is mg^p = 0.35 GeV. Left (right) panel corresponds to the = ^qcd ~ 0.2 GeV 
{'fng,v = 0.7 GeV) case {rug^^ is the gluon mass in the vacuum). 
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Figure 8; The difference between medium and the vacuum 0*^ order fractional energy loss for light 
(dotted curve), charm (dashed curve) and bottom quark (dot-dashed curve) is shown as a function 
of the thickness of the medium. The initial jet energy is E = 15 GeV. The curves are computed 
by assuming running coupling given by Eq. with different couplings in the medium and in the 
vacuum. Light quark mass is M = 0 GeV, charm mass M = 1.5 GeV, and bottom mass M = 
4.5 GeV. The gluon mass in the medium is rug^p = 0.35 GeV. Left (right) panel corresponds to the 
mg^u = Aqcd ~ 0.2 {mg^y = 0.7) GeV case {rug^y is the gluon mass in the vacuum). 
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Figure 9: The difference between medium and the vacuum 0*^ order energy loss for charm quarks 
is shown as a function of coordinate of particle production. The curves are computed by assuming 
running coupling given by Eq. m, with different couplings in the medium and in the vacuum. 
The thickness of the medium is L = 5 fm. The charm mass is M = 1.5 GeV and the gluon mass 
in the medium is mg^p = 0.35 GeV. Left (right) panel corresponds to the = h-qcD ~ 0.2 
("^' 5 ,'!; = 0.7) GeV case (mg^^ is the gluon mass in the vacuum). 
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Figure 10: The difference between medium and the vacuum 0*^ order energy loss is shown as a 
function of the gluon mass in the vacuum (mg^y). The curves are computed by assuming running 
coupling given by Eq. m, with different couplings in the medium and in the vacuum. The dotted 
curve corresponds to the light quarks, dashed to the charm and dot-dashed curve to the bottom quark. 
The initial energy of the jet is 15 GeV, and thickness of the medium is L = 5 fm. We take that light, 
charm and bottom quark masses are M = 0, M = 1.5 and M = 4.5 GeV respectively. 
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between medium and the vacuum energy loss is in the range of about ±15%, which may have a sizable 
effect on the pion suppression. We therefore conclude that, in order to obtain consistent predictions 
for pion suppression data, 1) the gluon mass in the vacuum has to be more accurately estimated and 
2) the transition radiation has to be taken into account. 

5 Net radiative energy loss dependence on transition radiation 

In this section we will use the medium induced radiative energy loss given in m to study how 
the difference between net radiative medium ± and the vacuum (AE'ial) energy loss 

changes when the transition radiation effects are included. 

5.1 Comparison between light and heavy quark medium induced radiative energy 
loss 
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Figure 11: On the left panel the 1'^* order in opacity fractional energy loss as a function of initial jet 
energy is shown for heavy and light quark jets. Assumed thickness of the medium is L = 5 fm. On 
the right panel the 1®* order in opacity fractional energy loss for a 15 GeV jet is plotted versus the 
effective static thickness L. Plasma is characterized by rug^p = 0.35 GeV and A = 1 fm. Dotted curves 
corresponds to light quarks while (dashed) dot-dashed curves corresponds to charm (bottom). 

For the purpose of further comparison, in this subsection, we show the 1®* order medium induced 
radiative energy loss for light and heavy quarks. The left panel of Fig. shows fractional energy 
loss to 1®* order in opacity(AFid)) as a function of initial jet energy. As in the previous sections, 
we assume running coupling given by Eq. m with rUg = rug^p = 0.35 GeV. We take L = 5 fm and 
A = 1 fm for the plasma parameters. We see that for 5 GeV jet, the finite mass effect leads to a 50% 
(90%) reduction of the energy loss for charm (bottom) quarks. On the other hand, for a 20 GeV jet, 
we see that the finite mass effect has almost no effect on charm quarks while it reduces the bottom 
quark energy loss by 50%. 

On the right panel of Fig. ^Jwe fix the jet energy to 15 GeV, and look at the fractional energy loss 
as a function of thickness of the medium. Wee see that charm and light quark energy loss dependence 
is similar, while bottom quark remains significantly different and close to the linear Bethe-Heitler 
form. This behavior is expected having in mind the left panel in Fig. IllL There we see that, for a 
15 GeV jet, the finite mass effect does not have a large influence on charm quarks, while it still has a 
significant influence on heavy bottom quarks. 
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5.2 


The net radiative energy loss for light and heavy qnarks in QCD medium 



central and right panels correspond to = 0, 0.2 and 0.7 GeV cases respectively. Dotted cnrves 
corresponds to light quarks, while (dashed) dot-dashed curves corresponds to charm (bottom). 


In this subsection, we first concentrate how the net radiative energy loss depends on the initial jet 
energy for „ = 0, 0.2 and 0.7 GeV cases. Figure IT^ was obtained by combining Figs. El and Q with 
the left panel of Fig.ITTl We use it to compare the net radiative energy loss results for light, charm and 
bottom quarks. We see that, depending on the gluon mass in the vacuum, the transition radiation may 
either further enhance (for > mg^p) or kill the “dead-cone” effect (for mg^^ < m-g^p). Additionally, 
we see that the transition radiation may have a significant influence on the net radiative light parton 
energy loss. For example, for the mg^y = 0 GeV case, the light quark energy loss is smaller than the 
charm quark energy loss. Additionally, in this case, the energy loss for all three types of quarks shows 
a weak dependence on the initial jet energy. On the other hand, in the rng^y = 0.7 GeV case, the light 
quark energy loss is a steeply decreasing function of initial jet energy. Additionally, this energy loss is 
significantly larger than both the charm and bottom quark energy losses. 



the thickness of the medium is shown for heavy and light quark jets. Initial jet energy is 15 GeV. Left, 
central and right panels correspond to rug^y = 0, 0.2 and 0.7 GeV cases respectively. Dotted curves 
corresponds to light quarks, while (dashed) dot-dashed curves corresponds to charm (bottom). 
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Figure ITHl is obtained by combining Figs. 0 and |S1 with the right panel of Fig. We use it to 
compare the thickness dependence of the net radiative energy loss results for light, charm and bottom 
quarks. In the < Aqcd case, we see that for the light quarks, introduction of the transition 
radiation may lead to the cancellation of the medium induced radiative energy loss for L <3 fm. This 
result infers that light partons may leave the medium practically unquenched if traveling the distances 
smaller than 3 fm. Similar energy loss thickness dependence was already observed in d, and the 
Fig. IT^mav point how to solve the thickness dependence puzzle posed by M- 
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Figure 14: The net radiative fractional energy loss {AE = + AE^^^ — AE^^c) as a function of 

the gluon mass in the vacuum {mg^^) is shown for heavy and light quark jets. Dotted curve corresponds 
to light quarks, while (dashed) dot-dashed curve corresponds to charm (bottom). Initial jet energy is 
= 15 GeV, and thickness of the medium is L = 5 fm. 
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Finally in Fig. dl we fix the jet energy and thickness of the QCD medium, and compare the rug^y 
dependence of the net radiative energy loss results for light, charm and bottom quarks. While bottom 
quark net radiative energy loss is independent on mg^y, light quark shows significant dependence on 
the on nig^y, as expected from Fig. ^1 For example, we see that, for the nig^y < Aqcd tbe light 
quark energy loss is smaller than charm’s. On the other hand, in nig^y ~ 1 GeV case, the difference 
between light and heavy quark energy loss is enhanced from 10% (see 15 GeV point on Fig. ini) to 
approximately 40%. 

Unfortunately, at the moment we do not know what value of rug^y would, most accurately, reflect 
confinement effects in the vacuum. However, based on Figs. II211141 we see that the rug^y > rug^p 
enhances the differences between the light and heavy quark energy loss results, and therefore would 
lead to a significant difference between the suppressions of light and heavy parton observables. On the 
other hand, rug^y < mg^p lowers the differences between the light and heavy quark energy losses, and 
would correspondingly lower the differences between the suppressions of light and heavy observables. 
Based on the most recent experimental results EZl, which suggest similar suppression results for pions 
(light partons observable) and single electrons (heavy quark observable), we expect that mg^y < Aqcd 
is the most appropriate value to approximate confinement in the vacuum. 

6 Conclusions 

A finite size medium, with dimensions on the order of the diameter of the coiiided heavy ion, is created 
in URHIC. Due to that, jets experience a transition from medium to the vacuum, which results in 


15 




additional energy loss, called transition radiation. Since the Ter-Mikayelian was computed under the 
assumption of an infinite medium m, in this paper we addressed the finite size correction to this 
effect. 

In |1()) we obtained that, due to the Ter-Mikayelian effect, the medium energy loss for charm 
quarks is reduced by 30 % compared to the vacuum case. We have showed that the finite size 
correction reduces this effect from 30% to 15%. The obtained result is intuitively unexpected, since 
the well known QED transition radiation calculations iniEni give a positive difference between the 
medium and the vacuum energy loss. The discrepancy between QED and QCD case results from the 
fact that in QED experiments a particle is produced far outside the medium, and has to cross two 
vacuum/medium boundaries in order to reach a detector, while in QCD experiment, the particle is 
produced inside the medium and crosses only one boundary. The QCD effect is therefore smaller by 
(approximately) the energy loss corresponding to one boundary crossing. 

Previously, there was a contradiction caused by the fact that the energy loss is infinite in the vacuum 
and finite in the medium, leading to the inhnite discontinuity between medium and the vacuum energy 
loss. This problem was long avoided by assuming the same zero mass for light partons in both medium 
and the vacuum. We here showed that this infinite discontinuity is naturally regulated by including 
transition radiation. To our best knowledge, the work presented here is the first consistent solution 
to this problem. 

Further, we showed, that for rug^^ < Aqcd, the light parton may not loose energy when traveling 
distances smaller than 3 fm. This result is similar to the one experimentally observed in [14| . Con¬ 
sequently, one of the future goals is to understand the impact of the transition radiation to the light 
observable’s suppression results. 

We note that our computations were done under the assumption of static medium of finite size L. 
Therefore, one of the interesting future problems is to study how these results are modified under the 
influence of the dynamically expanding QCD medium with continuously changing density at the edge. 
We expect that for a more gradual density change between the medium and vacuum, the difference 
between medium and vacuum 0*^ order radiation is reduced. 

Finally, we also note that our calculations considered only the radiative energy loss and did not 
take elastic energy loss into account. Recent computations by [28| show that the elastic energy loss in 
the QCD medium is negligible, which supports non-inclusion of the elastic energy loss into account. 
However, in m it was obtained that the elastic energy loss is significant. Consistent inclusion of the 
elastic energy loss into our radiative energy loss formalism is the subject of our future work. 
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A Gluon wave function in finite size QCD medium 

In this appendix we will derive the gluon wave function in a finite size QCD medium. As already 
stated in Section |2J in the spinless case the wave function of the emitted gluon with momentum k can 


16 


be written as (see Eq. ©) 


\{x) = e^,{k) ^>g(x) c, (22) 

where e(k) = [o, 2 |^,e] is the transverse polarization and c is the color factor of the radiated gluon. In 
finite size QCD medium, gluon mass (which is proportional to the temperature) becomes position 
dependent, and ‘hg(x) is the wave function that satisfies the Klein-Gordon equation with position 
dependent gluon mass 

In this section we will derive the wave function ^g{x) in the light cone coordinate system. To 
compute this wave function, we have to solve the fallowing equation 

(□ + mg{i&f)^g{x) = 0, (23) 

where in light cone coordinate system □ = Ad^+d^- —d^. Note that we have an extra factor of 4 which 
is the consequence of the coordinate transformations. Additionally, note that Klein-Gordon equation 
now becomes first order in x"*" and x“ ■ [ISlISIj. 

By assuming that the 4>g(x) has the fallowing form ‘hg(x) = (/>i(x+)()) 2 (x“) it becomes easy 
to show that 

</,2(x-) = e-'5^^-- (24) 


where k'^ is a constant. 

Klein-Gordon equation then reduces to 


dln(/)i(x+) .k^-|-mg(x+)^ 

dx+ * 2/c+ 




where k (x"'') = iL+lhsllLJ-, 


We can now easily obtain the solution for ())i(x'’') 


</>i(x+) 


X+ 

-f / dik-{i) 


(25) 


(26) 


which, together with Eq. (EH), leads to the solution of the gluon wave function in a finite size QCD 
medium 





[/c+x -\- j d^k (^)]+2k-x 
0 


(27) 


This solution is valid for arbitrary gluon mass functional dependence mg(x'*“), and is not limited 
to the static medium case which we consider in Section IH In a static medium, where mg{y.'^) is given 
by Eq. 0, the Eq. (EH reduces to 


where 


$^(x) = e{L 


^ ) I -kv )L -ik^x _ j^\ 

2 ’ ^ 2 ’ 




kv = [k'^ 1 


k+ 
k^ -|- m 


2 

'g,'v 


k+ 


k] 


is the gluon momentum in the medium and the vacuum respectively. 


(28) 


( 29 ) 
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B Computation of Mrad in the case when the jet is produced inside 
the medium 

In this appendix we will compute the amplitude of the diagram Mrad in the case when the jet is 
produced inside a static QCD medium of size L. To do that we start with the Eq. o, i.e. 

M'^ad ^ j ^4^^ Am(xi - Xo) u^(xi) (xi)$i(xi). (30) 

Here <h(xi) = is the wave function of the final quark with (on-shell) momentum p and H^(xi) is 
the wave function of the emitted gluon. Vertex function u^(xi) is given by u^(xi) = g'(xi)( — 5^), 

where ^(xi) is the running coupling constant. 

After replacing Eqs. dH)-® in Eq. we obtain (note that since xi > xq, we keep only the first 
term in Eq. ®) 

Mrad = J J(xo) d^Xi J ^((^1 “ Xo) + ) 

9{^i)(d^^ - ^a(xi) c 

{p' + pYe^ik) - xo)0((xi - xo)+) c. (31) 

In the static medium we can replace d?g(x) by Eq. (EHl), and the Eq. m reduces to (note x = xi— xq) 


< (Sy /(){x+)«(L - V) 

+ j d^x {-i)gy YiP+>^^-p')^ -L)}c, (32) 

where we used J(p -|- A;) = f (Yxq J(xo)e*(P+^)^°. 

We will first compute h = f (~‘^)9p e^^M/«p-p')^ 0(x'^)0(L — ^). Note that in the light 
cone gauge d^x = l/2dx+dx“dx, leading to 

Ii = /’l/2dx+dx-dxe5(P+'=-P')+^” e^(P+^^-P')'^^e-*(P+^-P'>^d(x+)0(L-—) 

{27rY J 2 

r2L 

= -9p + A:+ - p'+)d(p k - p') / dx+ eYP+^p-p')~ 

do 

1 _ pi{p+kp-p')~L 

= -2i5rp (5(p+-h/c+-p'+)d(p-hk-p')-w—7-7^ (33) 

{p + kp- p') 

Without loss of generality, we can take for the initial quark the plane wave state in the x-plane and 
set p' = 0. Then p = —k, and {p + kp — p')~ = Xp, where Xp is given by Eq. ® . By using this, the 
Eq. (IHHl) finally reduces to 

1 _ pi-XpL 

h = -2igp d(p+ + k^ - p'~^)S{p k - p')- (34) 
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In the same way I 2 = ( 2 ^ f i~'^)9v “ ^) reduces to 


I 2 = —2igp e +k'^ — p^’'')5(p + k — p') 


oiXvL 


Xv 


= —2igp d{p~^ + k^ — p~^)6{p + k — p 




Xv 


(35) 


After we plug in the Eqs. (1511) and (jlH) into Eq. (inU and use (p' ■ e) = ^ {x = ^), the Mrad 
reduces to 


Mrad = -2i J dp'~^(fp' J{p + k) {p ■ e) (5(p+ + k~^ - p'+) (5(p + k - p') 


^ ffp 1 - e^^p^ ^ gy e^^p^ 


P 


/+ 


Xp 


-) c 


, N e k 

= —2iJ{p + k) - 

X 


P'^ Xv 
' gp I- e^^p^ gy e*^r 

■n'+ 


Xp 


p'+ Xv 


(36) 


which is the Eq. (|SI) given in section [21 


C Computation of Myad in the case when the jet is produced outside 
the medium 

In this appendix we will compute the Eq. m in the case when the jet is produced at the distance Iq 
from the medium, i.e. at Xg" = —21q < 0. Then, for x = xi — xg, (where xi is the gluon production 
point), mg{x^) can be written as 

mg(x'^) = nig^y 6'(x+)6'(/o - ^) + rug^p e{^ - lo)0{lo + L- ^) + rug^y 9{^ - (L + Iq)) (37) 
By using Eq. (1571) . ‘hp(x) (see Eq. (1571) 1 reduces to 

^^(x) = eix+)e{lo - ^) + 0(^ - Zo)0(/o + L-^) 

^^-i(k--k-)L^-ikr^ (38) 

With the use of Eq. (i55|). the Eq. m reduces to 


M™. = l^-!^Ap + k)i2p'.e) 

{ <-')«” e'"’*''--’’'’* «(x+)9(io - 

+ ^iik--k-)lo ^ f ^4^ (_^) g*(p+fcp-p')x 0(^^ _ +L-^) 

[liry J 2 2 

+ J d^x {-i)gy e^iP+kr-p')^ ^(^ _ (^ + ;^))| ^ _ ( 39 ) 

By applying the same procedure as in Appendix B, My ad finally reduces to 
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Mrad 



J{p + k) {p' ■ e) 6{p'^ + k~^ — p'~^) 5(p + k — p') 



2i / dp' 


J{p + k) {p' ■ e) (5(p^ + k~^ — p'^) (5(p + k — p') 



g^ _ g^(Xp-^“l“X't'^o) 


P'+ Xv 


Pv e^(xpL+xvio) 


2iJ{p + k) — ^\—-{—- ^)e^x.io(^i _ ^ixpL^ 
X p'+ IXv Xv Xp 


(40) 


which is the Eq. (O given in section 01 
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